subroutine read_scalars()

use storage
use mpi_storage

implicit none

logical :: file_exists

inquire(FILE='scalars.dat', EXIST=file_exists)
if(file_exists) then
   if(n_rank == 0) write(*,*) 'Reading namelist scalars from scalars.dat'
   open(unit=133,FILE='scalars.dat')
       read(133,nml=scalars)
   close(133)
endif

end subroutine read_scalars
